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Abstract 

The magnetic moment fi of the p-meson is calculated in the framework of QCD sum 
rules in external fields. Bare loop calculations (parton model) give: fipart = 2.0. The 
contribution of operators of dimension 6 reduces this value: n = 1.5 ± 0.3. 
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1 Introduction 



Investigation of the static properties of vector mesons provides an important informa- 
tion about strong interaction of hadrons. In particular, the vector dominance hypothesis 
(VDM) supposes that the interaction of real or virtual photon with hadrons proceeds 
in such a way that the photon first transforms into vector mesons p, uo, (p, which then 
undergo interaction with hadrons. In the consistent lagrangian formulation of VDM it is 
assumed [0] (for review, see 0) that p-mesons are Yang-Mills vector bosons. In the frame- 
work of this hypothesis the p- meson magnetic moment is equal to 2 (in units eh/{2'mpc)), 
at least if strong interaction is neglected. 

The goal of this paper is to calculate the p- meson magnetic moment in QCD, using 
the method of QCD sum rules in external fields 0,0]. 

In paper the p-meson formfactors were found at intermediate momentum transfer 
by QCD sum rules. By extrapolation of the p-meson magnetic formfactor to the point 

= (outside the applicability domain of the technique) it was found that the p-meson 
magnetic moment p is close to 2. However, this result can not be considered as conclusive; 
the direct calculation of p in QCD in model independent way is still absent. The p-meson 
magnetic moment was calculated in the Dyson-Schwinger equation based models 0,0 
and in the framework of relativistic quantum mechanics [Q. 

Here we work in the limit of zero quark masses, a^-corrections are neglected. 

2 Phenomenological part of the sum rule 

We consider the correlator of two vector currents in the external electromagnetic 

field: 

n^.(p) = 1 1 d'xe'^^{T{3,{x)3tmF . (1) 

Here subscript F denotes the presence of the external electromagnetic field with strength 
Fp\ and is the vector current with p-meson quantum numbers: = U'~f^d. Its matrix 
element is 

{p'-\3,\0) = iml/gp)e,, (2) 

where rup is the p-meson mass, gp is the p-7 coupling constant, (7^/(47r) = 1.27, and is 
the p-meson polarization vector. 

In the limit of weak external field we consider only linear in Fpx terms in the correlator 
(1): 

H^, = Hj, + iVI^U^^^^F^^ . (3) 

We find magnetic moment from sum rule for the invariant function H(p^) at certain 
kinematical structure of ^pux^r (3). To obtain this sum rule, we calculate H at < as 
the operator product expansion series. On the other hand, we saturate dispersion relation 
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for n by the contributions of physical states. After equating of these representations the 
required sum rule appears. 

Therefore, first of all one should choose kinematical structure. 

The electromagnetic vertex of the p- meson has the following general form P]: 



(p(p + g,e'")|j^'|p(P,e")) = -e^'e^f {i2p + q)^gp^ - {p + q)pg^^ - p,gpJF^{-q') + 



- ^x-9p)^2(-g') + ^(p + q)pPa{'2p + g)x^3(-g') ) • (4) 

nip 



■el 
'X 

and Fi , F2 , F3 are electric, magnetic and quadrupole formfactors correspondingly, 



In (4) j'^ = CuU'y^u + edd'j^d is electromagnetic current, e„ , are u- and (i-quark charges 



Fi(0) = l, /i = l + F2(0), (5) 

fi is the p-meson magnetic moment. 

Using (2) and (4), we obtain for the (0|j^j|p) (p|j^|p) (p|j,y|0)e^ transition: 

-^E(o|j.l/)(/b?l/)(Plj.|o) 



ex 



r,r 



= ^^Y.K^^'^^yy^xiij.'^P + q)x9po - (P + (l)p3x'y - Pa9px)Fii-q'^) + (6) 

+ {9xpQa - gxaqp)F2{-q^) + ^(p + q)pPa{'2p + g)x^3(-g^)) ■ 

•nip J 

Here e^— photon polarization, r, r' are the p-meson polarization indices. Let us consider 
in this expression linear in q^ terms. We sum over p-meson polarizations, retain the 
antisymmetric over x, a part, introduce F^^- = i^e^qa — ^aQx) obtain for (6): 
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^4 . 11 

■^F^aUF2 + -Fi) — (py{p^gp„ - Pag^,x) - Pf^{Px9ua - PaQvx) 



2g2 2 

+ ^Fi^[p^{p^gp„ - Pag^^x) + Pp{Px9ua - PaQux)) + (^2 + Fi){gp^g^„ - g^agux)^ ■ 

Formfactor F3 does not give linear in q„ contribution. 
Thus, we choose the structure 

Pu{Px9pa - Pa9,^x) - Pt^{Px9ya - PaQux) ■ (7) 

In comparison with another possible structure, g^ix9i'(^~9i^(^9i^x^ i^) contains two additional 
powers of momentum in the numerator, which result in better convergence of the operator 
expansion series. 
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It should be noted here that, as follows from the vector current conservation, the 
antisymmetric over field indices x, a structure in H^u^a (3) is antisymmetric over p-meson 
indices fi, v too. 

From (5) one can see that ^2(0) + l/2Fi(0) = /i - 1/2. 

Thus, one should calculate the invariant function n(p^) at the structure (7) in H^uxrr- 
In the dispersion relation for 11 we use the simplest model of physical spectrum, which 
contains the lowest resonance and continuum. Phenomenological representation of 11 has 
the form: 

1 / 1 



Here dots mean the contributions of non-diagonal transitions (for example, 

ij.\p*) {p*\jx^x\p) {p\j'^\^) y where p* is the excited state with the same quantum numbers 
as p), function / represent continuum contribution and Sp is the continuum threshold for 
the p-meson. 

Retaining only the terms, which do not vanish after Borel transformation, we obtain: 



where C appears due to non-diagonal transitions. 

3 Calculation of the vector current correlator 

Now let us calculate n(p^), basing on the operator product expansion in QCD. 
The quark propagator in the external electromagnetic field F^^, in the fixed-point gauge 
XpA^ = , A^ = —l/2F^yXi, can be found in 0: 



327r2x2 



Here is the quark charge, a, {3 are spinor indices, a, 6 - color indices and (see 0) 
{qcTpxq)F = V '^T^CiCqxFpx {qq) , x is the quark condensate magnetic susceptibility. 

The expression for the quark propagator in the external electromagnetic and soft gluon 
fields has the following form in momentum representation |]^: 

Sfg = ^ (7A7r7p7aP - 2pA7r7<77p - 2pr7A7p7f7- 

-8ppPr9Xa . , ^ . ^ .\ 
P + ^gXrlalpP - ^9XT9paP) ■ 
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Here G"^ is the gluon field strength and are the color matrices. 
The contribution of the loop diagrams to 11 (p^) is equal to 



IGtt^ J (s — p^Y 

According to the quark-hadron duality, the continuum contribution in the interval of 
= —p^ from Sp to infinity is determined by the bare loop in this interval. Therefore, 
function / in (8) is constant: / = — 3/(167r^). 

The loop diagrams correspond to the operator of the lowest dimension Fpx. Operators 
of dimension 4 are absent. As was shown in 0, operator q(Dp'ji, — D^'jfj)q has opposite 
with respect to the electromagnetic field C-parity and can not be induced by them, while 
operator ep^^pxqj^'ypDxq vanishes due to equation of motion for massless quarks. 

There are a number of vacuum expectation values of operators of dimension 6: 
{qcrpxq)F{qq), {Gl^Gl^)Fp, and 

g{q[{GlxD. - "D^Glx) - {GlxD, - pG:x))lxt''q) f , 

e.upxgm'l^Dx + SAG;\)7e75t"g)F , (10) 

jikl ( /^i /^k /^l /^i /^k /^l \\ 

" y^fiX^Xp^pu ~ ^uX^Xp^pp.))F , 

where Dp is the covariant derivative and d*'^' are SU(3) structure constants. 

The diagrams, corresponding to the operator {G^^G^^)Fp^, have infrared divergence. 
We introduce the cut-off over transversal momenta A and obtain their contribution into 
n(p2): 

^]_(^r^\(^ ^ + ^"1 nn 

36 ^vr ^V2AV 6AV ' 
This divergence is probably cancelled by the contribution of field induced vacuum ex- 
pectation values (10) (for example, this was shown in [Q for dimension 4 operators and 
symmetric tensor field). Usually such vacuum expectation values can be calculated by 
constructing corresponding sum rule. But for operators (10) this approach is inapphcable 
because of their high dimension. 

The dominating contribution appears from no loop diagrams with hard gluon ex- 
change. In our case such diagrams contain the operator {q<Jpxq) F{qq) ■ They give: 

'^ 9^{qq?x ..^x 

9 / ■ 

It should be noted here that the quark condensate magnetic susceptibility x is negative. 

Collecting the expressions (9), (11), (12), one can find the operator product expansion 
part of the sum rule: 

n p =-7^/7 iW + o 6 ^ H^mT^ ~ + ~j ■ ^2 

loTT^ J [s — p^)"^ 9 p° 36 TT ^2X^p^ byp^ p°' 

The uncertainty in the contribution of the vacuum expectation values (10) is included 
into the errors. 



ds , , 

(9) 
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4 Results and discussion 



After Borel transformation 

( p2\n+l / W \ " 

B{M^)= lim ^-^-^(--^) , P^ = V>0 

we equate the phenomenological (8) and operator product expansion (13) parts of sum 
rule and obtain: 

f.-l- + CM^ = (l - e-^^/'^e-VM^- 



ml \ 



9M2 36 ^ TT 



M2 



1 

3A2 



M2) 



(14) 



C appears due to nondiagonal transitions. 

We use the following values of parameters: 
rrip = 0.77 GeV— the p-meson mass, 
(^2/(47?) = 1.27— the p-7 coupling constant, 
Sp = l.SGey^— the continuum threshold for p-meson, 
((a,/7r)G2) = 0.009 ± 0.007 Ge\/^-the gluon condensate 
g^{qqy = (0.28 ± 0.09) x 10"^ GeF^-the quark condensate 
X = —(5.7 ± 0.6) GeV^'^—the quark condensate magnetic susceptibility 

= 0.8 GeV'^—the cut-off over transversal momenta. 

First of all, let us consider the contribution of the bare loop (and continuum). It is 
given by the first term in (14). In |jl2| the following relation for Qp can be found: 



m 



A'K'^m?p ^ 



s,/NPyml/M' 



e p' 



(15) 



Substituting (15) into (14) and omitting for a while the terms with quark and gluon 
condensates, one can obtain very simple answer: 
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We see that in the parton model approximation the p-meson magnetic moment is equal 
to 2. This result agrees with prediction of the vector dominance hypothesis. 

Now let us analyze the whole equation (14). In order to find the value of magnetic 
moment, we approximate the right-hand side of (14) (see fig.l) by a straight line in the 
interval 0.9 Ge^^ < < 1.3 GeV^ and find it ordinate at zero Borel mass. 

Thus we obtain: 

p=1.5. (16) 

The contribution of the operators of dimension 6 to this value does not exceed 20%. 
The contribution of the terms, which contain A^, is not more than 20% of the total 
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contribution of dimension 6 operators. That is why variation of within the interval 
0.6 GeV^ < < 1.0 GeV"^ does not change the value of magnetic moment. 

The variations of the values of the quark and gluon condensates within the given limits 
change the value of magnetic moment by < 10% each. 

The uncertainty in the value of the quark condensate magnetic susceptibility results 
in the error about few percent in the value of the magnetic moment. Variation of the 
continuum threshold for the p-meson Sp in the reasonable limits gives the same effect. 

Supposing that the contribution of vacuum expectation values (10) does not exceed 
50% of that from diagrams with hard gluon exchange (12), we obtain after collecting all 
uncertainties: 

/i= 1.5 ±0.3. 

This is our final result. 

In Ijl^ it was shown that the approximation procedure is correct (nonlinear terms can 
be safely neglected), when /i > CM^. In our case CM^//i ?a 0.2 0.3. 

Thus we find that in parton model (bare loop) approximation the p-meson magnetic 
moment ^part = 2, whereas the nonperturbative interactions decrease this quantity by 
a quarter: /i = 1.5 ± 0.3. It is important to mention that all accounted operator prod- 
uct expansion corrections are negative, i.e. results in decrease of p in comparison with 
fJ'part = 2. Since the effective values of the Borel parameter are about IGeV'^, one 
may expect that perturbative corrections are remarkable and can reach ~ 20%. 

The value of the p-meson magnetic moment was calculated in a number of papers 
within the Dyson- Schwinger equation based models. In 0] the value /i = 2.69 was found. 
In 1^ several results are compared, and the values of /x lie between 2.5 and 3.0. Relativistic 
quantum mechanics model gives [§] A* = 2.23 ± 0.13. Unfortunately, while a^-corrections 
are not calculated (we plan to do it in the next paper), it is hard to say with certainty if 
this discrepancy is real or not. 

The author is grateful to B.L. loffe for posing the problem and valuable discussions 
and to A.G. Oganesian for helpful discussions. 

The work is supported in part by grants CRDF RP2-2247, INTAS 2000 Project 587 
and RFFI 00-02-17808. 
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Figure 1: The right-hand side of equation (14) R{M'^) as the function of M^. 
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